Quantum phase transition induced by Dzyaloshinskii-Moriya in the kagome 

ant iferromagnet 
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We argue that the S = 1/2 kagome antiferromagnet undergoes a quantum phase transition when 
the Dzyaloshinskii-Moriya coupling is increased. For D < D c the system is in a moment-free phase 
and for D > D c the system develops antiferromagnetic long-range order. The quantum critical 
point is found to be D c ~ 0.1J using exact diagonalizations and finite-size scaling. This suggests 
that the kagome compound ZnCu3(OH)eCl3 may be in a quantum critical region controlled by this 
fixed point. 



In the search for materials realising a spin-liquid 
ground state, one has to face the presence of small 
anisotropic interactions of spin-orbit origin. Such in- 
teractions that break the full rotation symmetry of the 
Heisenberg model, reduce the quantum fluctuations and 
may tend to induce magnetic phases at low tempera- 
tures. The recently discovered spin 1/2 copper oxide 
ZnCu3(OH)6Cl;?i that has the geometry of a kagome 
lattice may be a good candidate for a spin-liquid. 2 Ex- 
perimentally no apparent freezing of the magnetic mo- 
ments has been found down to very low temperatures 
despite strong Heisenberg interactions. Exact diago- 
nalizations of the Heisenberg model predict indeed a 
non-magnetic state with no magnetic moment ~> 7 How- 
ever, smaller interactions of spin-orbit origin are cer- 
tainly present. In particular those of Dzyaloshinskii- 
Moriya symmetry^ are expected when the magnetic 
bonds have no inversion center, which is the case of 
ZnCu3(OH)6Cl3.— An immediate question is to what ex- 
tent they affect the non-magnetic phase. 

Dzyaloshinskii-Moriya interactions have been first 
invoked^ in the context of ZnCu3(OH)eCl3 to ex- 
plain the enhancement of the spin susceptibility at low 
temperatures. 5 NMR measurements of the local suscepti- 
bility have provided a different interpretation in terms of 
the presence of defects in the structure^ a result consis- 
tent with a direct fit of the susceptibility^ and corrob- 
orated by theoretical calculations. 12 A more direct evi- 
dence of anisotropy has been obtained by paramagnetic 
resonance^ A Dzyaloshinskii-Moriya coupling of order 
0.08 J was needed to explain the linewidth, 13 which is typ- 
ical of cuprates. It is therefore a relatively small correc- 
tion but may be of crucial importance in highly frustrated 
systems. In fact it is known that an infinitesimally small 
Dzyaloshinskii-Moriya interaction in the classical kagome 
favors long-range Neel order with a Q = propagation 
vector and 120° orientation of the spins. 14 Spin- wave cor- 
rections renormalize down the magnetic moment but do 
not suppress it for D = O.U. 15 So it was unclear how 
this Neel phase could be reconciled with experimental 
observations. 

In this letter, we show that the proper inclusion of 
quantum fluctuations leads to a phase transition from a 



Neel state to a moment-free phase at a quantum critical 
point that we estimate to be D c ~ O.U. This therefore 
resolves the contradiction for ZnCu 3 (OH) 6 Cl3 and, fur- 
thermore, suggests that its magnetic properties may be 
strongly influenced by the proximity of the quantum crit- 
ical point. In particular power-law scalings were observed 
in the dynamical susceptibility- and NMR relaxation 
tirtieS) 10 i 16 and have been interpreted so far in terms of 
critical spin liquid states for the kagomei 17 ' 18 The present 
results indicate that the power-laws of ZnCu3(OH) 6 Cl3 
may well originate in the present fixed point. In order to 
identify the phases in an unbiased way and locate D c , we 
have performed exact diagonalizations of small clusters. 
At D = O.U, there is a clear emerging low-energy tower 
of states that collapses onto the ground state like 1/N, 
thus signaling a broken-symmetry phase. We have calcu- 
lated the Neel order-parameters as function of D/J and 
have found that they vanish for D < O.U, thus leading 
to a moment-free phase and D c ~ O.U. 




FIG. 1: (color online). The kagome lattice with the allowed 
Dzyaloshinskii-Moriya Interactions. The orientation of the 
bonds specify the order of the operators in Si x Sj. m is a 
mirror plane. 

The model we are considering is based on the symme- 
tries of the two-dimensional kagome lattice (Fig. [T])^ 
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where nn stands for nearest neighbors and S, is a S = 1/2 
quantum spin on site i. Thanks to the lack of inver- 
sion symmetry at the middle of each bond, interactions 
of Dzyaloshinskii-Moriya symmetry are allowed between 
nearest neighbors. According to Moriya's rules^ there 
are components of perpendicular to the planes, of 
strength D, staggered from up triangles to down trian- 
gles; and in-plane components that point towards the 
center of each triangle of strength D' (see Fig.QJjiii There 
is a priori no simple relation between D and D' . It is 
true that if the Cu-(OH)-Cu plane was a mirror plane 
of the crystal structure, then the D vector would be 
perpendicular to it. In fact the proton of (OH) breaks 
that symmetry— and since the perturbation is expected 
to be strong it is difficult to relate D and D' . For 
ZnCu3(OH)6Cl3, J has been estimated from the suscepti- 
bility to be about 170-190 K. 9,1 - Interplane couplings are 
thought to be smaller because copper ions are far away 
in the c direction and separated by sheets of zinc and will 
be neglected in the following. According to electron spin 
resonance measurements, D = 15 K (D ~ 0.08J) and 
D' = 2 K (£>' ~ 0.01.7)3 

We now map the model ([T|) onto a simpler model that 
restores a U(l) symmetry up to terms of order D' 2 /J. 
The in-plane components vectors sum up to zero when 
going around a triangle and are therefore reducible to a 
D' 2 1 J term (for small D')^- by appropriate rotations of 
the spin operators^ We shall neglect these second-order 
terms not only because they are smaller but also because 
there are other (symmetric) exchange anisotropies at the 
same order that we have not included. Therefore in the 
rotated frame, the model has only the original D com- 
ponent along z with the same strength and has a U(l) 
rotation symmetry about this axis: it is this symme- 
try that, as we shall show, is going to be spontaneously 
broken. In the following we shall present the results of 
the numerical diagonalization of (HJ for systems of size 
N = 21,24,27,30,36 (for the A = 24 to 27, we have 
considered 2 different cluster shapes). The dimension of 
the largest Hilbert space is ~ 7 x 10 8 . We have started 
with a Dzyaloshinskii-Moriya coupling strength fixed to 
D = 0.1J (and J= 1). 

Symmetry breaking in the thermodynamic limit. We 
have calculated several low-energy levels of (TTJ) in each 
sector of the total magnetization, denoted by S z , for dif- 
ferent cluster sizes. For A — 36, the energy levels are 
shown in Fig. [2j The spectrum is qualitatively different 
from the exact spectra obtained at D = (see refs. [?})■ 
Here we clearly see a band of low energy levels (red line in 
Fig. [J) well separated from higher energy states. It forms 
a so-called tower of states which energy is very well de- 
scribed by a quadratic term S 2 /2Im (Fig. [2]). The slope 
1/(2Jjv) is fitted for all available A and shown in the inset 
of Fig. [21 It clearly extrapolates to zero in the thermody- 
namic limit like 1/A. In this limit one can then form a 
superposition of these eigenstates with different S z , thus 
forming a macroscopic state with a preferred in plane di- 
rection. This therefore shows that the system breaks the 




rl i I i I i I i I 

5 10 15 20 



FIG. 2: (color online). Low-energy levels of the Heisen- 
berg kagome lattice with Dzyaloshinskii-Moriya interactions 
for N = 36. The inset gives the slope of the lowest energy 
states versus S%, denoted by 1/(27jv), as a function of 1/N. 
(75 = 0.1.7). 



rotation symmetry in the thermodynamic limit. In this 
standpoint, the degeneracy of the tower of states is then 
a natural consequence of the invariance of the Hamilto- 
nian in rotations about the internal axis defined by D: a 
given macrosopic state pointing in the u direction in the 
transverse plane is degenerate with all its transforms in 
any rotation about D. In addition, the wave- vector of the 
lowest energy state in each S z sector is Q = for all clus- 
ters, so that we can safely conclude that the system will 
not break translation invariance in the thermodynamic 
limit. Furthermore, because of the continuous broken- 
symmetry, we expect a long wavelength Goldstone mode 
with energy varying like k oc A -1 / 2 in two dimensions. 
Unfortunately, the first allowed k ^ wave-vectors are 
not small enough on these clusters to observe the indi- 
vidual long wavelength states directly in the spectra. In 
fact, indirect confirmation of a Goldstone mode will be 
given below in the scalings of the energy of the ground- 
state and correlations. From the spectra, we can now 
extract the uniform susceptibility at T=0 by using the 
expression of the energy, S z /2In — HS Z . The suscep- 
tibility of the ground state In is proportional to A (as 
shown in Fig. [2]) and the susceptibility per site x — In /A 
is found to be x = 0.144 ± 0.002 for D/J = 0.1 in the 
thermodynamic limit. 

Neel order-parameters. To test for the possibility to 
have a sublattice magnetization (or Neel order) in the 
thermodynamic limit, we have calculated the spin-spin 
correlations in the (finite-size) ground state. We distin- 
guish between correlations in a plane perpendicular to 
D, (0\Sf a S* b + S% a S] b \0), and the correlations along D, 
(0\S* a SJ b \0) . As expected from the easy- plane character 
of the Dzyaloshinskii-Moriya interaction, the latter re- 
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mains much smaller and especially at large separation. 
To study the in-plane ordering, we define the Fourier 
transform of the spin-spin correlations: 



^(Q) = ^E eiQ ' (R< " 



R, 



(0\Sf a S%\0) (2) 



This is a 3 x 3 matrix that turns out to be peaked at 
Q = 0. The largest eigenvalue at Q = corresponds 
simply to the 120° in-plane orientation of the magnetic 
moments within the unit-cell. The prefactor 24 is cho- 
sen so that the largest eigenvalue is f in the perfect 
Neel state. For a quantum magnet with a finite Neel 
order-parameter, the eigenvalue should extrapolate to a 
finite value smaller than f (because of quantum fluctu- 
ations), with a well-defined finite-size scaling. The lat- 
ter can be predicted from the existence of a low energy 
spin wave with wavevector varying like N~ x / 2 in two di- 
mensions, and is of the form N^ 1 / 2 for the correlations 
and iV~ 3 / 2 for the energy per sitei^i Fig. [3] shows the 
largest eigenvalue versus N- 1 ' 2 for different D/J. First 
we see that the scaling is obeyed both for the correlations 
(Fig. [3]) and for the groundstate energy per site (inset 
of Fig. [3]) . Second, the intercept for infinite system-size 
gives the order-parameter, noted m 2 AF . Starting from 
large D/J (squares) and reducing D/J we see that the 
order-parameter decreases from m 2 AF ~ 0.326 to zero at 
D c ~ 0.1J (circles). For D < D c , it is no longer pos- 
sible to find a finite order-parameter: the negative ex- 
trapolated values reflect the breakdown of the scaling at 
some length-scale and short-range correlations develop 
instead. In Fig. 01 we summarize the behavior of the Neel 
order-parameter tuaf, as a function of D/J. The order- 
parameter decreases continuously to zero, so that the 
transition is compatible with a second-order phase tran- 
sition. It seems difficult though to extract an accurate 
value for the critical exponent, given the error bars re- 
sulting from finite-size scaling. The present data are com- 
patible with the mean-field behavior, m 2 AF oc (D — D c ) 
but the exponent could also be smaller. In any case, 
D c ~ 0.1J clearly appears as the critical point separat- 
ing the Neel phase (Q — 0, 120° in-plane orientation of 
the spins from a phase with no static moment. It is a 
more accurate estimate than that of spin-waves (that is 
three times smaller)^ because here all quantum fluctu- 
ations have been taken into account. 

In order to have a simple understanding for the occur- 
rence of a Neel phase in the phase diagram, we discuss the 
limit of large D/J, for which simple physical arguments 
can be used. For this we start with a trimerized version 
of the kagome latticed On a triangle, three Heisenberg 
quantum spins 1/2 form a quartet and two degenerate 
low-energy doublets. The Dzyaloshinskii-Moriya interac- 
tion lifts the degeneracy and selects the doublet with the 
(vector)-chirality opposite to D. We now assume that D 
is large and ignore the higher doublet (and quartet). We 
define a pseudo-spin o; L — 1/2 for the lowest doublet on 
each triangle. The inter-triangle Dzyaloshinskii-Moriya 
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FIG. 3: (color online). Finite-size extrapolation of the largest 
eigenvalue of S a b(Q = 0), showing (i) a scaling in N' 1 ^ 2 ; 
(ii) the decrease of the sublattice moment from m\ F = 0.326 
(D = 1, J = 0) to zero (D = 0.1J) and negative values 
(D — 0.06; 0.08J). Inset: finite-size scaling of the ground 
state energy in N~ 3/2 for D = 0.1 J. 
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FIG. 4: The Neel order-parameter as a function of D/J, 
obtained from extrapolations to infinite size. The system has 
no magnetic moment for D < D c ~ 0.1J and has Neel order 
for D > D c . 



interaction (that is now supposed to be large) gives an ef- 
fective interaction for the pseudo-spins that we can write 



H ' = -q\ D \ E 



■aVa^ + eij ■ (pi x a 3 ) (3) 



<»j> 



where i,j are now sites on a triangular lattice, ey is a 
unit-vector along +z when turning anticlockwise around 
an up triangle of the triangular lattice. is essen- 
tially a ferromagnetic in-plane interaction with a size- 
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able effective Dzyaloshinskii-Moriya interaction. The lat- 
ter gives magnetic frustration and the problem of solv- 
ing ([3]) is in principle as complicated as the original 
problem. For three triangles, the ground state of ([3]) 
is the Q = doublet state (of pseudo-spins), made of 
(I HI) + I TIT) + I ITT) )/V3 and its time reversal coun- 
terpart. Similarly, the classical ground state of ([3]) is the 
Q = in-plane ferromagnetic state with e.g. (<rf) = 1/2 
on all triangles. In terms of the original spins, this cor- 
responds to long-range Neel order in (Si a ), of the Q = 0, 
120° form. The reduction of the in-plane moment can 
also be estimated. For the three-triangle ferromagnetic 
doublet given above, one can form a superposition of the 
two states that gives a moment along x, for instance. The 
moment is then found to be a fraction vtiaf — 4/9 ~ 0.44 
of the full moment. This gives a simple explanation for 
the Q = Neel phase found by exact diagonalization 
in the limit of large D/J. Since we know from previous 
works that the D — phase is non-magnetic*^ we natu- 
rally expect at least one quantum critical point between 
the two phases. 

We conclude that the S — 1/2 kagome antiferromag- 
net has a quantum critical point at D c ~ 0.1J separating 
a moment- free phase (D < D c ) from a Neel phase with 
a sublattice moment (D > D c ) (or a weak net moment 
if D' 7^ 0, see [HI). This is clearly compatible with the 
absence of a static moment in ZnCu3(OH)6Ci3j^£ be- 
cause the coupling extracted from ESR, D = 0.08 J, 13 
is smaller than the critical coupling^ Furthermore this 
estimation places this compound very close to the quan- 



tum critical point. This raises the issue of the origin of 
the power-law behaviors observed experimentally when 
the temperature is decreased ^ 10 ' 16 They have been in- 
terpreted so far in terms of critical spin liquids ; 17 ' 18 or 
free 5 = 1/2 impurities^ Here we suggest a different 
intrinsic interpretation in terms of the proximity with 
the present critical point. The prediction of the criti- 
cal behavior and low-temperature scalings is not an easy 
task, however. In the view of the present uncertainties 
about the D = phase and the other possibilities of va- 
lence bond crystal states^ or intermediate phases such 
as a spin ncmatic state, for instance, it is indeed diffi- 
cult to ascertain what the effective low-energy quantum 
field theory is. It is interesting to note though that the 
same ordered phase was found as an instability of the 
field theory describing the algebraic spin liquid^ 7 - but 
here we have given a finite critical value for D c . At finite 
temperatures such a critical point will open a quantum 
critical region that may have consequences on the mag- 
netic properties of ZnCu3(OH) 6 Cl3. It is also of course 
of particular interest if the compound could be driven 
across the transition by applying an external pressure or 
magnetic field, for instance. 
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